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ABSTRACT 

The cosmic microwave background (CMB) polarisation and the 21 cm line fluctuations 
are powerful probes of cosmological reionisation. We study how the cross-correlation 
between the CMB polarisation (E- modes) and the 21 cm line fluctuations can be used 
to gain further understanding of the reionisation history, within the framework of 
inhomogencous reionisation. Since the -E-mode polarisation reflects the amplitude of 
the quadrupole component of the CMB temperature fluctuations, the angular power 
spectrum of the cross-correlation exhibits oscillations at all multipoles. The first peak 
of the power spectrum appears at the scale corresponding to the quadrupole at the 
redshift that is probed by the 21 cm line fluctuations. The peak reaches its maximum 
value in redshift when the average ionisation fraction of the universe is about half. On 
the other hand, on small scales, there is a damping that depends on the duration of 
reionisation. Thus, the cross-correlation between the CMB polarisation and the 21 cm 
line fluctuations has the potential to constrain accurately the epoch and the duration 
of reionisation. 
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1 INTRODUCTION 



The history of the cosmological reionisation is one of the open problems of modern cosmology. Questions like what causes the 
reionisation and ho w it proceeds are intima tely related to the evolution of matter density fluctuations and to the formation of 
the first structures ( Barkana fc Loebl 2001 ). The lack of observational data makes it quite difficult to ans wer these questions 
today . Up to now, the available probes of the reionisation epoch are just a few, e.g. the Gunn-Pet erson test rtGunn fc Peterson 
1965) and the cosmic microwave background (CMB) radiation, in particular its polarisation |Zaldarriagalll997h . The ob- 



servation of the Gunn- Peterson effect, which probes the amount of neutral hydrogen, suggests that full reionisation was 
accomplished by z ~ 6 ( Fan et al. 20061 ). The CMB polarisation, in turn, measures the optical depth of the Thomson scat- 
tering from the l ast scattering surface to the present epoch. The WMAP 3 years data provide an optical depth r ~ 0.09 
1 Page et al.ll2007l ) suggesting that reionisation began around z = 10. 



In addition to those two probes, the observati on of fluctuations of the 21 cm line background is expected to be one of 
the most promising methods to study reionisation 1 Madau et al. 1997 ; Tozzi et al. 2000l ; Ciardi fc Madau 20031 ). The 21 cm 
line corresponds to the energy of the hyperfine spin flip transition of neutral hydrogen atoms. After decoupling of the baryon 
temperature from the CMB temperature around z ~ 300, neutral hydrogen atoms absorb or emit a 21 cm line, depending 
on their temperature. Since the 21 cm line is redshifted by the cosmological expansion, we can obtain redshift slices of the 
Universe by choosing the frequency of observation. It is thus possible to follow the evolution of the i ntergalactic medium before 
and during reionisation directly. The 21 cm line observations have recently gained a lot of attention IjLoeb fc Zaldarri aga 2004; 
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Furlanetto et al hood I Tas hiro fc Sugivam al l200d : ICooravll200d ) and several projects are designed for measuring the fluctuating 
line background (e.g. MWaQ, LOFARq SKaI). 

In the near future, we will therefore have different direct and complementary probes of reionisation. We thus expect that 
cross-correlations between them will provide even more information than their respective auto-correlations. This is the case 
especially f or the cross-correlat i on between the CMB temper ature anisotropies and the 21 cm line fluctuations on large scales 
(£ ~ 100 



>r the cross-correlation between the L-MrJ temperature anisotropies and the 21 cm line fluctuations on large scales 
Alvarez et"al"]l20od : lAdshead fc Furlanettdl2008T ) and on small scales (I > 1000 - ISalvaterra et alJliood : ICoorav 



20041: Isiosar et al.ll2007h . Since the observations of the 21 cm lines are expected to be dominated by foregrounds and noise 



jjelic et alj 20081 ) . calculating the cross-correlation is indeed a good way to extract information on reionisation. In particular 
on large scales, the 21 cm line fluctuations are corr elated with the CMB Doppler temperature anisotropies which are due to 
the motion of ionised baryons. Alvarez et al. ( 2006h showed that the maximum amplitude of the cross-correlation is reached 
at t he redshift when the ionised fraction is one half, which potentially provides a way to constrain double reionisation models 
fe.g. lCer]|2003h . 

In the present study, we examine the cross-correlation between the CMB i?-mode polarisation (a more effective probe 
than the temperature anisotropies) and the fluctuations of the 21 cm line background on large scales. Measuring _B-mode 
polarisation is one of the present challenging observations for CMB experiments. However, projects like Planck and CMB-pol 
will provide precise CMB polarisation maps in the near future. Similarly, cosmological observation of 21 cm line background is 
still belongs to the future, but LOFAR and MWA are being completed and will soon provide us with maps of neutral hydrogen 
at high redshifts. Hence, it is useful to investigate what information on reionisation we can extract from the cross-correlation 
of those maps. Indeed, we show here that the cross-correlations are a potentially powerful tool to constrain the history of 
cosmic reionisation. This article is organised as follows. In Sec. II, we provide the analytic formula of the cross-correlation 
between the CMB _E-mode polarisation and the fluctuations of the 21 cm line background. In Sec. Ill, we show and discuss 
the results of the cross-correlation. Sec. VI is devoted to the conclusions. Throughout the paper, we use WMAP 3-year values 
for the cosmologica l parameters, i.e. h = 0.73 (Ho = h x lOOkm/s/Mpc), T = 2.725K, h 2 Q b = 0.0223 and h 2 Q m = 0.128 
( Spergel et al.ll2007l ) for a flat cosmology. We also set the speed of light to c = 1. 



2 CROSS-CORRELATION BETWEEN CMB £-MODE POLARISATION AND 21 CM LINES 
2.1 The CMB £-mode polarisation 

The CMB polarisation is best described by the Stokes parameters, Q and U. Sin ce the Stokes para meters have spin 2, they 
can be decomposed onto the plane waves and spin-2 spherical harmonics ±2 Y™ l|Hu fc Whitelll997t ). 



(Q±iU)(r,,x,n) =J2(-i)'J ^1 j ( E t } ± exp(ifc ■ x) ±2 Y/™(n). (1) 

£ m 

Here, m = 0, ±1, ±2 correspond to the scalar, vector and tensor types. In this paper, we only consider m = modes, because 
they dominate other polarisation modes at the scales we are interested in. 

The CMB polarisation components, Ej, and B^ , are obtained from the Boltzmann equations for CMB. According to 
those, Bp =0 for all is and the evolution of £^ ' is given by the following integral solution, 



E?\w,k) _ 3 / (I + 2)! n° T m Mk(VQ-v)) (2) 



where P*- ' is the m = source term due to Thomson scattering, f is the differential optical depth for Thomson scattering in 
conformal time, given by f = n e a^a with the electron number density n e , the cross section of Thomson scattering <jt and the 
scale factor a normalised to the present epoch, and r is the optical depth between the conformal times 77 and 770- The source 
term P^ ' is related to the initial gravitational potential <E>o via the transfer function Ds(fc, 77), 

P (0) =D E {k, v )$ . (3) 

We provide the detailed calculation of the transfer function in the appendix. 
We also define the _E-modes in the direction n on the sky by 



1 http:/ /web. haystack. mit.edu/array/MWA 

2 http://www.lofar.org 

3 http://www.skatelescope.org 
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2.2 21 cm line fluctuations 



The o bserved brightness temperature of the 21 cm lines in a direction n and at a frequency v is given as in iMadau et al 
(1993) by 

T-21 



T 2l {h;u) = 



- (T s — Tcmb ) (Vobs ,n(rj — r] obs ) ) • 



(5) 



(1 + Zobs) 

Here, the conformal time 77 b s refers to the redshift z obs that satisfies v = 1/21 /(l + z obs ) with ^21 being the frequency 
corresponding to the 21 cm wavelength. The optical depth for the 21 cm line absorption is 



T21 ~ 8.6 x 10~ d xh- 



Ts 



n b h 2 
0.02 



0.15 



10 



(6) 



where xh is the fraction of neutral hydrogen, and T s is the spin temperature given by the ratio of the number density of 
hydrogen atoms in the excited state to that of hydrogens in the ground state. 

Equations ([5J and ([6]) show that the observed brightness temperature of the 21 cm lines will reflect the baryon density 
fluctuations 5 b and the fluctuations of the ionised fraction <5 X . We can rewrite Eq. ((5} in the linear approximation for both 8 b 
and 8 X as 



Tax (ft; 1/) = [l-x e (l + S x )](l + 5 b )T , 
where we x e — 1 — xn is the ionised fraction and 
0.15 \ / 1 + Zobs 



To = 23 



fl b h 2 
0.02 



On/! 2 



10 



1/2 



mK. 



(7) 



(8) 



Here, we assume that T B is much larger than the CMB temperature which is valid soon after the beginning of the reionisation 
l|Ciardi fc Madaull2003h . 

After performing the Fourier transform of Eq. (J7J), we obtain the 21 cm line fluctuations, 



ST 2 i{h-v) 



T 2 x(ft;z/)-To 
dk 3 



To 



t [(1 - x e )S b — x e 5x] exp [—ik ■ ((rj — ri obs )n)]. 



(2tt) 3 

Therefore, the 21 cm line fluctuation map at the frequency v can be described by 



STo 



L (ft; v) = T ^j ^3V4tt(2^ + 1) [(1 - x e )(l + F/i 2 )5 b - x e S x ] j e (k( Vo - r, obs ))Y e ° (n) . 



where we applied Rayleigh's formula to Eq. (|9}, 
exp(-ifc • x) = J2 VMM + l)(-ifj e (kr)Y e °(h), 



(9) 



(10) 



(11) 



where n is measured in the coordinate system where fc is e$. In Eq. (j 10 p - we included the factor (1 + F/i 2 ) to account for 
the enhancement of the fluctuation amplitude due to the redshift distortion (Kaiser effect (|Kaiserlll987l )) on the 21 cm line 
fluctu ations where [s, = k ■ n and F — d\n g / d\n a with g(a) the linear growth factor of baryon fluctuations (jBharadwai fc Ali 
2004 ). 



2.3 Cross-correlation between 21 cm and i5-mode fluctuations 

The angular power spectrum is defined as the average of the spherical harmonic coefficients ai m over the (21 + 1) m- values, 



n _ (Km| 2 ) 

^-L. (27+T)' 

where the ae m are defined for an arbitrary sky map f(n) as 



(12) 



(13) 



From Eqs. ((4]) and (110|l . we can thus derive the angular power spectrum of the cross-correlation between i?-modes and 21 cm 
line fluctuations 



Cf 



4ttT / dri 



dk 3 (EUv, k) [(1 - Xe)(l + Fn 2 )8 b - x e S x ] ) 



(27T) 



(2£ + l) 



-(1 — x e )P$s b — X e P$s : 



jdHvo - Vv)) 

Mk(r)o - Vobs)) ji(k(r]o - ??)) 



(k(r) - v)) 2 



(14) 
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where P<i>s b and P<s>s x are the power spectra of the cross-correlation between the gravitational potential and the baryon density 
fluctuations, and between the gravitational potential and the fluctuations of the ionised fraction, respectively. In the second 
line, we use F(fi 2 ) = 1/3 for the matter dominated epoch. 

The cosmological linea r perturbation theory pr ovides the relation between the baryon density fluctuations and the initial 
gravitational potential (e.g. Kodama fc Sasaki 1984 ). 5b(k,rf) — k 2 D b (k, if)$o(k) where Db is the transfer function of baryons. 
The power spectrum P<s>s b can thus be written in terms of the initial power spectrum of the gravitational potential P$ as 



P$s b = k D b (k,r])P<s>. 



(15) 



Although the reionisation process is not well-known, we can reasonably expect that there is a relation between the fluctuations 
of the ionised fraction and that of the matter density. Ionising sources are formed in dense regions. They ionise the surrounding 
medium with an efficiency that depends on the density of the medium. This translates into two possible cases: One where 
ionis ed fluctuations and matter over -densities coin c ide, a nd one where ionised fluctuations and matter density are anti-biaised 



(e.g. Benson et al. 2001). Following Alvarez et al. ( 20061 ). we assume that the fluctuations of the i onised rate are associated 



with the matter density contrast using the Press-Schechter description (|Press fc Schechterl 1 19741 ) . As a result, the power 
spectrum P<s s x is given by 



-xn In xn[bh 



f]D m {k, V )k 2 P 9> 



(16) 



where D m is the transfer function of matter (both dark and baryonic), &h is the average bias of dark matter halos more 
massive than M m i n , the minimum mass of the source of ionising photons 

h = l + J- - f , r , (17) 

where a(M) is the variance of the density fluctuations smoothed with a top-hat filter of the scale corresponding to a mass M, 
and /coil is the fraction of matter collapsed into halos with M > Mmin- In this paper, we choose M m i n such that the halo virial 
temperature is 2Vir(Afmin) = 10 4 K. This choice corresponds to the assumption that the ionising sources form in dark matter 
halos where the gas cools efficiently via atomic cooling. The parameter / describes the reionisation regime we are interested 
in. For / = 0, we are in the "photon-counting limit" case where recombinations are not important and where the progress 
of the reionisation depends on the number of ionising photons only. The over-dense regions contain more collapsed objects 
which are sources of ionising photons. Therefore, in this case, ionisation in over-dense regions is easier than in under-dense 
regions. On the contrary, / = 1 indicates the "Stromgren limit" case where ionisation is balanced by recombination. Although 
the over-dense regions contain more sources of ionising photons, the recombination rate in over-dense regions is higher than 
in under-d ense regions. Hence, over-dense regions in the / = 1 case have a lower ionised fraction than in the / = case (for 



details, see lAlvarez et al.ll200' 



Finally, in order to calculate Eq. (I14p . we need the evolution of the mean neutral hydrogen fraction for which we use 
a simple parameterisation based on two key quantities, the reionisation epoch (defined as the redshift at which the ionised 
fraction equals 0.5), z la , and the duration of reionisation, Az, 

1 



x H (z) 



1 + exp[— (z — z IC )/Az] ' 
In the following we will vary both z T , 



(18) 



and Az, and investigate their effect on the cross-correlation signal. 



3 RESULTS AND DISCUSSIONS 

As one can see in Eq. (|14p . the angular power spectrum Cf~ 21 of the cross-correlation is a function of several key parameters 
that encode the prominent details of the reionisation history. Consequently, and as we show here, the global shape of the 
Cf~ 21 as well as the amplitude and redshift dependence of its first peak are features that can be traced back to the reionisation 
characteristics. 

First it is worth noting that because of the double integration of the spherical Bessel functions in Eq. (|14|l , the angular 
power spectrum reflects the source term P(k obB ), namely the quadrupole term of the CMB temperature anisotropy, at z obB 
where fc b a satisfies fc b 3 = l/{rio ~ Vobs). Accordingly, the angular power spectrum of the cross-correlation exhibits its first 
peak at a multipole £ < 10 which corresponds to the angular separation of the quadrupole at z obs . The oscillations at higher 
is (> 10) are due to the free streaming of the quadrupole at redshifts higher than z obs . 

We first investigate the effect of the duration of reionisation Az. We present in the left panel of Fig. [1] three cases: One 
corresponding to a instantaneous reionisation Az — 0.01 and two with longer durations (Az = 1 and 2). We compute the 
angular spectra of the cross-correlation between 21 cm line fluctuations and the CMB P-mode fluctuations for those different 
durations but for the same given epoch of reionisation, z rc — 10, an observing redshift, z b s = 10, and in the case of the 
"photon-counting limit" (/ = 0). The difference in the duration of reionisation induces both a minor and a major effects 
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z„ = 10, Az = 0.01 
z* = 10, A- = 1 
z». = 10, A; = 2 
E-mode polarization 



^ = 10, As = 2 - 
4. = 10, A 2 = l - 
a, = 10, Az = 0.01 ■ 




Figure 1. The left panel shows the angular power spectra of the cross-correlation between .E-mode polarisation and 21 cm line fluctu- 
ations. We set z rc = 10, z b s = 10 and / = for all curves. The dotted line represents the angular spectrum for Az = 0.01. The solid 
and dashed lines are for Az = 1 and Az = 2, respectively. For comparison the thin solid line, shows the angular power spectrum of the 
CMB E-modes for z ro = 10 and Az = 1. The right panel is the evolution of the visibility function fe~ T in the same reionisation model 
as the left panel. The dotted, solid and dashed lines are for Az = 0.01, Az = 1 and Az = 2, respectively. 




Figure 2. The angular power spectra of the cross-correlation between .E-mode polarisation and 21 cm line fluctuations. The dotted line 
is for z = 12. The solid and dashed lines are for z = 10 and z = 8, respectively. We use Az = 1, z Q b s = 10 and / = for all curves. 



on the spectra. The minor effect is the observed shift of the peak position to higher Is with increasing durations. Since a 
longer duration increases the contribution of the visibility function re _T in Eq. (|14[) . which we plot in the right panel of 
Fig. [1] the contributions from redshifts larger than z b s accumulate in the calculation of the integration. However, spherical 
Bessel functions have a steep cutoff after r; b s so that the contributions from epochs earlier than r/ b s are negligible. Therefore, 
additional effective contributions come from only the former epoch rj < rjobs> and make the peak shift toward higher £ modes. 

The major effect of Az exhibited in the plot (left panel of Fig. [TJ is the damping at high Is with increasing durations of 
reionisation. This occurs again due to the increase of the amplitude of the visibility function at z > z Te . At high is, the source 
term P' ' exhibits rapid oscillations. Since the tail of the visibility function is larger, the range of the integral is larger, and it 
thus encompasses more oscillations. As a result, many more cancellations occur which damp the overall power. We note that 
the damping depends on the duration of the reionisation which is an advantage of observing the cross-correlation signal. As a 
matter of fact, it is difficult to see differences among CMB polarisation spectra for different durations, because the Thomson 
optical depth is not strongly sensitive to this parameter. However the plots in the left panel of Fig. [T] clearly show that the 
damping of angular spectrum of the cross-correlation is sensitive to the duration of reionisation making it a potential tool for 
constraining this parameter of the reionisation history. 

We now vary the reionisation redshift z IC keeping the other parameters fixed, and setting z b s = 10, Az = 1 and / = 0. 
We find that increasing z IC shifts the position of the first peak to small scales (Fig. [2| . The reason is the same as the one 
invoked for the shifts observed in the left panel of Fig. [I] When reionisation happens early, the integration of the .E-mode 
polarisation in Eq. (|14[) takes into account more contribution from higher redshifts. This induces a shift of the peak position 
to higher £ modes. 



Figure|3]shows the evolution of the first peak amplitude as a function of z bs, the redshift of the 21 cm line emi s sion, for the 
three above-mentioned durations of reionisation Az. Here, we set z TC = 10 and / = 0. As shown by Alvarez et al. ( 2006h using 
Eq. (118[) , the amplitude of the cross-correlation of the CMB Doppler temperature anisotropies and the 21 cm line fluctuations 
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Figure 3. The amplitude of the first peak of the cross-correlation signal as a function of redshift z b s - The dotted, solid and dashed 
lines indicate the amplitudes for Az = 0.01, Az = 1 and Az = 2, respectively. For all plots, we set z re = 10 and / = 0. 





Figure 4. The evolution of the first peak amplitude for different z rc . The dotted, solid and dashed lines represent the evolution for 
z r c = 12, z r c = 10 and z ro = 8, respectively. In this figure, we assume Az = 1 and / = 1. 



reaches its maximum value at z IC = 10 where the ionised fraction is about half. The same is true in the case of the cross- 
correlation of 21 cm line fluctuations and _E-mode polarisation. The reason for this is that the CMB polarisation is produced 
in the ionised regions whereas 21 cm line fluctuations are associated with neutral regions. We note that the signal in Fig. [3] 
is smaller at lower redshifts than at higher redshifts. This is due to the fact that the 21 cm signal vanishes once reionisation 
is completed, i.e. at redshifts smaller than z la . For higher redshifts (z b s > z rs ), the low scattered quadrupole component 
of the E- modes is suppressed like in the lower redshifts part, but the integration in the E- modes picks up contributions 
from the reionisation epoch and the iJ-modes are correlated with a larger contribution of the 21 cm signal. As a result, the 
amplitude of the cross-correlation does not decline as steeply as for z a ^ s < z TC . The difference between the curves is easily 
interpreted. In the Az = 2 case, full reionisation takes longer to complete than in the Az = 1 case. Therefore, neutral regions 
persist longer, broadening the skewed bell-shape of the curve centered on ;z b s = z rc . The slight difference in the peak value 
is due to the cumulative effect of the integral in Eq. flD - This is particularly well illustrated by the case of "instantaneous 
reionisation" (Az — 0.01). It is worth noting that our assumption, T s 3> T cm b, is formally not valid before the beginning of 
reionisation. In order to evaluate the cross-correlation signal correctly, we should rather take into account the full evolution 
of T s with redshift in the calculation of the 21 cm line fluctuations. However, by the time the ionisation fraction has risen 
abo ve the per cent level, efficient Compton and X-ray heating will guarantee that the spin temperature be well above T cm b 
(e.g. IPritchard fc Furlanettoll2007r i. 



The evolution of the first peak amplitude for different z rc is shown in Fig. [4] As expected, the amplitude reaches its 
maximum value at each z TC . Since the 21 cm optical depth is higher at high z b s , the amplitude of the 21 cm line fluctuations 
at higher redshifts is higher than that at lower redshifts. Therefore, and as exhibited by the comparison of the amplitudes for 
the different z le , the maximum value of the first peak cross-correlation amplitude is larger in the early reionisation case than 
in the late reionisation case. This again is the same evolution that shows up in the po wer spectrum of the cross-correlation 
between the CMB temperature anisotropy and the 21 cm line fluctuations (cf. Fig. 2 of Alvarez et al. 2006[ ). The evolution, 



with the observed redshift, of the shape and the amplitude of the first peak of cross-correlation between the 21 cm line 
fluctuations and Fv-mode fluctuations should allow us to constrain both the duration and redshift of reionisation. 
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redshift redshift 

Figure 5. The left panel shows the evolution of the first peak amplitude for two extreme values of the parameter /. The solid and 
dashed lines are for / = and / = 1, respectively. We set z rc = 10 and Az = 1 in all plots. The right panel shows the evolution of the 
first peak amplitude in the case of a double reionisation model. The solid line (left axis) is for the first peak amplitude at given redshift. 
The dotted line (right axis) indicates the ionisation fraction in the double reionisation model. 



Ultimately, the parameter / is not independent of the epoch z re of reionisation and its duration Az. However, making that 
dependency explicit would require the very detailed description of the reionisation processes that we are still lacking today. 
We can nevertheless illustrate how the cross-correlation power spectrum changes for different global reionisation scenarios. 
The left panel in Fig. [5] shows the variation of the evolution of the first peak amplitude with the parameter /. There is no 
difference of the peak position in both spectra for the "photon counting limit"/ = and the "Stromgren limit" / = 1 but 
the amplitudes are different. When / increases, recombinations become more and more important. Accordingly, over-dense 
regions in the case of / = have higher ionisation fraction than in the case of / = 1 as mentioned in Sec. 12.31 Since the 
fluctuations of the neutral fraction are defined by Sn — —S x , the absolute value of the 21 cm line fluctuations for / = is 
higher than for / = 1. Therefore, the absolute value of the amplitude of the CMB iJ-mode anisotropics and the 21 cm line 
fluctuations correlation for / = is higher than for / = 1 as shown in the left panel of Fig. [5] 

Finally, we also calculate the cross-correlation between the 21 cm line fluctuations and the _E-mode CMB anisotropies in 
the case of a double reionisation. For illustrative purposes, we chose a simple parametrisation of the mean ionised fraction 



xAz) = 1 : ; — — + A exp 

V ; l+cxp[-(z-z rc )/Az] y 



(z — Zfl ra t) 2 



(19) 



which is the simple reionisation case formula with z TC — 8 and Az = 1 on top of which we added a Gaussian with A = 0.8, 
Zftrst = 12, and Az| rst = 2 that encodes a first passage through an eighty per cent reionised universe at a redshift of twelve. 
These parameters have been adjusted so as to reproduce the same value of the optical depth (r = 0.08) as the single reionisation 
model shown in the left panel of Fig. [5] The redshift evolution of x in this toy model is represented by the dashed line in 
the right panel of Fig. [5] (read right axis for values of x e ). In the same figure, the solid line shows the corresponding evolution 
of the first peak of the cross-correlation (left axis). Every time the ionised fraction takes the value 0.5 there is a peak in the 
cross-correlation signal. In the case of our model for double reionisation, we observe three peaks, i.e. the ionisation fraction 
crosses the 0.5 line thrice. However, the two peaks at z b s ~ 8 and 11 are very close to each other. If two consecutive crossings 
of the 0.5 line occur very soon after each other, the c orresponding peaks in the cross-correlation signal will be essentially 
indistinguishable. Similarly to what has been found by Alvarez et al. 1 20061 ). this behaviour is a very interesting property as 
it implies that we can use the "trajectory" of the evolution of the peak amplitude to constrain the reionisation model. 



4 CONCLUSION 

We have examined the cross-correlation between the -E-mode CMB polarisation and the fluctuations of the 21 cm line 
background. In particular, we calculated the angular power spectrum of the cross-correlation on large scales using a simple 
parametrisation of reionisation, focusing explicitly on the redshift z la of its occurrence and on its duration Az. We also 
considered the case of double reionisation for which we obtained specific signatures. 

The angular power spectrum of the cross-correlation traces the quadrupole component of the CMB temperature anisotropies 
at the redshift where the 21 cm line is observed. Both the amplitude and the position of the first peak depend on the reion- 
isation epoch, as well as on the redshift z bs probed by the 21 cm observations. We showed that the signal peaks when z b s 



fluctuations of the 21 cm line background (see Alvarez et al. 2006), that is when the 21 cm observations probe the epoch 



at which the ionised fraction becomes one half. When probing higher redshifts, the polarisation and 21 cm cross-correlation 
signal decreases slowly, whereas it falls off rapidly when lower redshifts z \, B are considered. 
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The location of the first peak in the (z bs - Cf~ 21 ) plane is essentially independent of Az. Therefore, tracking the 
evolution of the first peak amplitude by scanning the sky in z b s bins would allow to determine z re . Moreover, the width of the 
bell-shaped curve describing the evolution with z D b s of the Cf~ 21 first peak depends on Az. The duration of the reionisation 
is thus also accessible through this observable quantity. 

We also showed that, in addition to the amplitude and the position of the first peak, the duration of reionisation induces 
a significant damping of the oscillations on small scales. This strong dependence on the duration shows specifically in the 
cross-correlation power spectrum as it does not arise in the auto-correlations power spectra. 

These significant features, which are the oscillations and the maximum at the redshift when the ionised fraction are one 
half, is explained by Eq. (|14|l which is an exact analytic expression. They are robust and independent of the details of the 
model. Assumptions are only introduced to relate the gravitational potential to the baryon density contrast and to the ionised 
fraction in Eq. (|14[1 . For the — 5& relation, we assume linear perturbation theory which is valid on large scales (I > 100) and 
at high redshift (z ~ 10). For the (j)-5 x relation, we need to assume a reionisation model. We choose a simple model based 
on th e halo model. In o ur redshift range when the reionised fraction is around one half, as shown in many numerical studies 
(e.g., Iliev et al. I 20061 )). the typical size of ionised bubbles is smaller than the instrumental resolution. In this case, we can 



as sume that the ionis ed fraction can be set by the number of ionising sources (for a detailed discussion see the appendix 
in Alvarez et al. 20061 ) . Numerical simulations that provide detailed descriptions of the reionisation history can test those 



assumptions. 

In the present study, we assumed the observation of the 21 cm line ideal. In practice, it amounts to consider that the 
frequency window function is the delta function 5(v — ^obs), implying that observations are done at the specified frequency 
l>obs only. Of course, for real observations, the finite instrumental frequency width has to be taken into account. For too broad 
window functions, this is susceptible to induce an additional damping of the higher I oscillations. We thus checked this effect 
in the case of the LOFAR interferometer that has a frequency width of 4 MHz. As a result, we did not obtain any additional 
significant damping of the oscillations for I < 200. 

Both CMB and 21 cm signals suffer from the foreground contaminations. Some of them will be correlated like the 
synchrotron emission from our Galaxy. Recent s t udies have shown that the cleaning of foreground were under control in both 
CMB and 21 cm observations ( Bock et al. 20061 : Jelic et al. 20081 ) . In addition, the multi frequency information on the CMB 



will allow the cross-correlation signal to be even less sensitive to the correlated foregrounds. 

Finally, due to the cosmic variance that limits CMB measurements on large scales, it may actually be difficult to pin-point 
the first peak with great accuracy. Moreover, ground-based telescopes designed for observations of the 21 cm line background 
cannot cover the whole sky. Those two facts are not important limitations for our study, since the damping of oscillations in 
the cross-correlation power spectrum on intermediate scales (10 < I < 100), which are well enough sampled, depends strongly 
on the reionisation evolution. Actually, even if the reionisation duration is so long that we cannot detect the oscillations due 
to the strong damping, it is still possible to obtain a lower limit for the duration. We will study this in more detail with 
numerical simulations in the future. 
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APPENDIX A: THE CMB £-MODE POLARISATION 



The source term of i?-mode polarisation in Eq. ([5J (|Hu fc Whitelll997f ). 

p (0) ( V , k) = e<°> („, k) - V6Ei 0) („, k), (ai) 

where is the quadrupole moment of the CMB temperature anisotropies. 

The CMB temperature anisotropies are evaluated by the following integral solutions obtained from the Boltzmann 
equations, 

ef\vo,k) p> 



^= fv-^ 4 o) (v)jf\k( V0 -„)), 

1 Jo „ 



21 + 

where the jjp are expressed in terms of the spherical Bessel functions je{x), 



Je 0) ( x ) = k{x), j^V) = j'i(x), j [ e Z> (x) = ^{3j" (x) + jt(x 



(2), 



and are source terms written as 



S, 



(o) 



(o) 



Co) 



tvb + fc*\ 



s<°> = - 

2 10 



(A2) 



(A3) 



(A4) 



where Vb is the velocity of the baryon fluid, 3> and 9 are the scalar-type metric perturbations in the conformal Newtonian 
gauge. The cosmological line ar perturbation theory p rovides the evolution of $ and * and the initial condition of Oo in the 
matter dominated epoch as (IKodama fc Sas aki 1984) 



$ : 



_9_ 

10 



*(0), 



O (0) = -*(0), 
5 



(A5) 



where *(0) is the initial metric perturbation. 

Before recombination, photons and baryons are tightl y coupled by Compton scattering. Hence, the tight-coupling ap- 
proximation is a strong tool for solving Eqs. ((2J) and (|A2[> . iHu fc Sugivamal |l995j) provided a useful analytic expression of 
the monopole component in this approximation, 



e (»7) + $(v) = [©0 (0) + $(0)] coskr a (rj) 



k_ 

71 Jo 



<V[$(t/) - tf (j/)] sin[fcr s (77) - kr s (r)% 



where r s is the integration of the sound speed c s = 1/ \/ 3(1 + R), 
1 



dr\- 



V3(l + R) 

The dipole and quadrupole components, and the baryon velocity are written as 



9i 



-(00 + *), e 2 = -^--9i, 



VB = ©1 



k K ~" ' ~" ~' 9 f" 
Accordingly, the quadrupole component of Z?-modes is given by 

V6, 



-02 



(A6) 



(A7) 



(A8) 



(A9) 



After recombination, the tight-coupling approximation is no longer valid. Using Eqs. (|A5|l - (|A9|) as initial conditions, we 
solve Eqs. ([2j and (|A2[) numerically in order to obtain the transfer function DE(k,rj) in Eq. 



